EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



CERN-PH-EP /2007-029 
DAPNIA-07-150 

9 July 2007 

Higgs boson searches in 
CP-conserving and CP-violating 
MSSM scenarios with the DELPHI 

detector 

DELPHI Collaboration 



Abstract 



This paper presents the final interpretation of the results from DELPHI on 
the searches for Higgs bosons in the Minimal Supersymmetric extension of the 
Standard Model (MSSM). A few representative scenarios are considered, that 
include CP conservation and explicit CP violation in the Higgs sector. The 
experimental results encompass the searches for neutral Higgs bosons at LEPl 
and LEP2 in final states as expected in the MSSM, as well as LEP2 searches 
for charged Higgs bosons and for neutral Higgs bosons decaying into hadrons 
independent of the quark fiavour. The data reveal no significant excess with 
respect to background expectations. The results are translated into excluded 
regions of the parameter space in the various scenarios. In the CP-conserving 
case, these lead to limits on the masses of the lightest scalar and pseudoscalar 
Higgs bosons, h and A, and on tan /?. The dependence of these limits on the top 
quark mass is discussed. Allowing for CP violation reduces the experimental 
sensitivity to Higgs bosons. It is shown that this effect depends strongly on the 
values of the parameters responsible for CP violation in the Higgs sector. 
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1 Introduction 

This paper presents the final interpretation of the Higgs boson search results from 
DELPHI in the framework of representative scenarios of the Minimal Supersymmetric 
Standard Model (MSSM). With respect to the previous MSSM interpretation published 
in Ref. [1], this analysis uses an enlarged set of experimental results, updated calculations 
of MSSM radiative corrections and covers more scenarios, including models with CP 
violation in the Higgs sector. 

As compared with the Standard Model, the MSSM has an extended Higgs sector with 
two doublets of complex Higgs fields, leading to five physical Higgs bosons, of which three 
are neutral. 




Figure 1: Main production processes of MSSM neutral Higgs bosons at LEP. Left: 
associated production of a Z and a Higgs boson, which must be one of the CP-even 
scalars (h or H) if CP is conserved or any Higgs boson (Hi, H2, H3) in the contrary case. 
At LEPl, the intermediate Z is on-shell and the final Z is off-shell, while it is the reverse 
at LEP2. Right: pair-production of neutral Higgs bosons. If CP is conserved, one of 
them must be CP-even (h or H) and the other one is the CP-odd pseudo-scalar A. If 
CP is not conserved, the pair can be any couple of different scalars among Hi, H2 and 
H3. The intermediate Z is on-shell at LEPl. 

If CP is conserved, two of the three neutral Higgs bosons are CP-even. They are 
denoted h, for the lighter one, and H. The third one is a CP-odd pseudo-scalar, denoted 
A. In e"'"e~ collisions, the dominant production mechanisms are the s-channel processes 
described in Fig. 1, that is the associated production of a Z and a CP-even Higgs boson 
and the pair production of either CP-even boson together with the CP-odd scalar. These 
processes are complemented by additional t-channel diagrams in the final states where 
a CP-even Higgs boson is produced with neutrinos or electrons, which proceed through 
W"'"W~ and ZZ fusions, respectively. These diagrams and their interference with the 
HiZ process have an impact on the production cross-section at masses around the HjZ 
kinematic threshold. At LEP 2 energies, the only significant effect is from W^W^ fusion 
which doubles the neutrino HjZ cross-section at the kinematic threshold. Finally, charged 
Higgs bosons, H"*" and H~, are produced in pairs through a diagram similar to that in 
Fig. 1, right, via exchange of a Z boson or a photon. 

Although CP is conserved at tree level in the MSSM, radiative corrections can in- 
troduce CP violation through stop and sbottom loops, leading to changes in the neutral 
Higgs boson sector [2]. If CP is not conserved, the three neutral Higgs bosons are no 
longer pure CP eigenstates but mixtures of CP-even and CP-odd components. They are 
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usually denoted Hi, H2 and H3, in increasing mass. The main production mechanisms 
are the same as in the CP conserving case, except that, a priori, any scalar can be pro- 
duced in association with a Z boson or through W"'"W~ and ZZ fusions, and any couple of 
different Higgs bosons can be pair-produced. The main phenomenological difference with 
respect to the CP-conserving case lies in the strength of the couplings of the Z boson to 
the Higgs scalars. In significant regions of the parameter space, CP violation turns off 
the otherwise dominant coupling between the Z boson and the lightest Higgs boson. In 
that case, if none of the other processes of Fig. 1 are possible (due e.g. to kinematics), 
the dominant Higgs boson production mechanism at LEP becomes the Yukawa process 
of Fig. 2. Of the two phases of LEP, only LEPl has a significant sensitivity to this 
process. In the Standard Model, the corresponding cross-sections are negligible, e.g. a 
fraction of a pb for a few GeV/c^ Higgs boson. In the MSSM, these can be enhanced by 
up to three orders of magnitude with respect to their Standard Model values, leading to 
detectable signals which become valuable in the case of CP violation. 




Figure 2: Additional production process of MSSM neutral Higgs bosons at LEP. The 
radiation of a Higgs boson off a Z boson decay fermion gives a detectable signal only at 
LEPl. This signal is exploited in the case of CP violation. 

The decay properties of the Higgs bosons are moderately affected by CP violation, at 
least in the range of masses accessible at LEP, that is up to masses around 100 GeV/c^ [2]. 
In most of the MSSM parameter space of the scenarios studied hereafter, the three neutral 
Higgs bosons decay mainly into the pair of heaviest fermions kinematically permitted, 
even if CP is not conserved. Below the threshold, a Higgs boson would decay into 

77 or e^e~ pairs with a significant lifetime. Above the fi~^fi~ threshold, the lifetime is 
negligible and Higgs bosons decay at the primary vertex. Up to a mass of 3 GeV/c^ the 
main decays are into pairs and also into hadronic channels with a large proportion 

of two-prong final states. Above 3 GeV/c^ the dominant decays are successively into cc, 
r+r~ and finally bb pairs for Higgs boson masses above 12 GeV/c^. Besides these decays 
into fermions, there are also regions of the parameter space where one neutral Higgs 
boson can undergo cascade decays to a pair of Higgs bosons, as for example h —>■ AA 
if CP is conserved or H2 —>■ HiHi in the contrary case. In some cases, especially if 
CP is not conserved, this mode dominates over the decays into SM particles. In the 
scenarios considered in this paper, charged Higgs bosons have a mass above 60 GeV/c^ 
and decay either into the pair of heaviest fermions allowed by kinematics, that is into 
cs or rz/ pairs, or into a W* and a light Higgs boson, e.g. into a W*A pair if CP is 
conserved. Finally, these scenarios do not allow neutral or charged Higgs boson decays 
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into supersymmetric particles such as sfermions, charginos or invisible neutralinos. Note 
that searches for neutral Higgs bosons decaying into invisible products were performed 
at LEP, as reported in Ref. [3]. 

The different decay channels define the topologies that were searched for to cover the 
MSSM parameter region kinematically accessible at LEP energies. These topologies 
are described in Section 2. Section 3 summarizes the definition and techniques related 
to confidence levels used in the statistical interpretation of the searches. The eight CP- 
conserving MSSM benchmark scenarios studied in this analysis are presented in Section 4 
and the results obtained in these scenarios when combining all searches are given in 
Section 5. Similarily, the CP-violating scenarios and the corresponding results are covered 
in Sections 6 and 7. The top quark mass has a significant impact on the properties of 
the Higgs bosons (e.g. mass spectrum of the neutral Higgs bosons, CP-violating effects). 
Results are thus derived for several values of this mass, namely: mtop = 169.2, 174.3, 179.4 
and 183.0 GeV/c^, which were defined after the measurement of the top quark mass at 
the Tevatron, run I [4] . Of the two values close to the present experimental measurement 
of 

"^top — 170.9 ± 1.1 ± 1.5 GeV/c^ [5], 174.3 GeV/c^ gives the most conservative results 
and thus was chosen as a reference in most of the exclusion plots and to quote absolute 
mass and tan/3 limits. Readers interested in similar analyses at LEP are referred to 
Ref. [6,7]. 

2 Search channels 

The different analyses performed to search for neutral and charged Higgs bosons in 
the whole LEPl and LEP2 DELPHI data samples are summarized in Table 1 which 
lists the final states, mass ranges, integrated luminosities and the references for more 
details about the selections and their performance. Two channels, the T"'"r~bb signal at 
LEPl and the (h — > AA — > cccc) (Z — > qq) signal at A masses below the bb threshold, 
were analysed for this paper, using selections already published. The efficiencies and 
the references for the selections can be found in the Appendices 1 and 2 of this paper. 
In the Tabic, the notations h and A which label the different analysis channels must 
be understood as generic notations for any pair of neutral Higgs bosons that could be 
produced in each of the production processes hsted in the Table. As an example, the 
hZ analyses, originally designed to search for the CP-even h boson in CP-conserving 
scenarios, can be applied to search for the second CP-even Higgs boson, H, as well as 
for any of the three Higgs scalars in CP-violating scenarios. It must be noted that the 
kinematic properties of the signal processes are only shghtly affected by CP-violation, 
since, when CP is not conserved, the production processes still proceed through the CP- 
even and CP-odd components of the neutral Higgs bosons, as explained in Ref. [7]. The 
same topological searches can thus be applied whether CP is conserved or not. 

As compared with our previous publication [1], the following changes were introduced 
in the experimental results used. The MSSM interpretation in Ref. [1] relied only on 
searches performed at LEP2 at masses above 12 GcV/c^ in Wh in the hZ process, with 
either direct or cascade decays, and above 40 GcV/c^ in mh, ruA in the hA channels, 
with only direct decays of the Higgs bosons. The corresponding channels have their ^/s 
values in bold characters in Table 1. Scans of the MSSM parameter space were thus 
restricted to rriA above 12 GeV/c^ and assumed the published LEPl limits-*^ to be valid. 
Including all LEPl results, which have a sensitivity starting from vanishing h and A 
masses, and the additional LEP2 searches of Ref. [21], whose sensitivity in the hA mode 

'^mh>44 (46) GeV/c^ when mh is above (below) the AA threshold [14] 
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bbbb, r+T-bb 


> 130. 


158.0 


2d 


[19] 


192-208 


T+r-bb 


> 120. 


452.4 


2d 


[20,1] 


192-208 


bbbb 


> 80. 


452.4 


2d 


[20,1] 


189-208 




> 8. 


570.9 


Id 


[21] 


189-208 


bbbb 


> 24. 


610.2 


no 


[21] 


189-208 


hadrons 


> 8. 


610.4 


mix 


[22] 



hZ or hA with h — > AA cascade 



91 Z ^ qq 

91 (AA ^ V°VO) (Z ^ any but r+r-) 

91 (AA ^ 77) (Z —* any or A ^ 77) 

91 (AA 4 prongs) (Z any or A ^ 2 prongs) 

91 (AA hadrons) (Z — > i/p or A — > hadrons) 

91 (AA T+T-T+T-) {Z^ VV or A ^ r+r-) 

161,172 (AA any)_(Z qq, vv or K ^ any) 

183 (AA ^ bbbb) (Z ^ qq) 

192-208 (AA ^ bbbb, bbcc, cccc) (Z qq) 

192-208 (AA cccc) (Z ^ qq) 

189-208 (AA ^ bbbb) (Z ^ qq or A ^ bb) 



< 0.21 


16.2 


no 


[15] 


< 0.21 


9.7 


no 


[15] 


< 0.21 


12.5 


no 


[15] 


> 0.21 


12.9 


no 


[15] 


> 0.21 


15.1 


no 


[15] 


> 3.5 


15.1 


no 


[15] 


> 20. 


20.0 


Id 


[17] 


> 12. 


54.0 


Id 


[18] 


> 12. 


452.4 


2d 


[20,1] 


> 4. 


452.4 


2d 


A.2 


> 12. 


610.2 


no 


[21] 



ffh or fTA Yukawa production 


91 


bb(h 


r+r-), bb(A ^ r+r-) 


4. - 50. 


79.4 


no 


[21] 


91 


bb(h 


^ bb), bb(A ^ bb) 


11.-50. 


79.4 


no 


[21] 


91 


r+r- 


(h T+T~), r+r-(A t+t~) 


4. - 50. 


79.4 


no 


[21] 


H+H- 


189-208 


cscs. 


CSTVr , W*ATVr , W*AW*A 


> 40. 


610.4 


2d 


[23] 


189-208 


T^V^ 


T~Vr 


> 40. 


570.8 


Id 


[23] 



Table 1: List of signals expected from MSSM Higgs bosons that were searched for in the 
DELPHI data sample. Indicated for each signal are the centre-of-mass energy, final state, 
analysed mass range, integrated luminosity, level of discriminant information inchided in 
the confidence level estimates (none, one- or two-dimensional) and the reference where 
details of the analysis are pubhshed. Here h and A denote any neutral Higgs boson allowed 
to be produced in each of the indicated production processes. The mass range applies 
to rrih for hZ production, to mh+mA for hA production, to rriA for h AA processes, 
to the Higgs boson mass for either Yukawa process and to for H+H^ production. 

When no upper bound is given, the limit imposed by kinematics or vanishing branching 
fractions must be understood. 
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complements that of the two other sets of results, allows scans of the MSSM parameter 
space to be performed with no restriction on masses. Moreover, some of the analyses 
of Ref. [21] cover production processes which are neghgible if CP is conserved but are 
enhanced by CP violation, such as Yukawa processes or the production of T~^T~T~^r~ 
final states. Adding the searches for neutral Higgs bosons decaying into hadrons of any 
flavour [22] is expected to provide sensitivity in scenarios where the Higgs boson decays 
into bb would vanish. As their mass coverage starts at low mass, these analyses also 
increase the experimental sensitivity to Higgs bosons below the bb threshold, a region 
otherwise covered only by analyses of subsets of the LEPl data. Finally, the charged 
Higgs boson searches [23] help in a few CP-conserving scenarios in the low tua region 
where the charged bosons are kinematically accessible at LEP2. 

Moreover, our previous interpretation was dealing only with the production of the two 
lightest Higgs bosons, the h and A scalars in CP-conserving scenarios. In this analysis, 
the production of the third boson, if kinematically accessible, is also accounted for, which 
can lead to a significant gain in sensitivity in restricted areas of the parameter space. 
In CP-conserving scenarios, this leads to including the HZ and HA signals besides the 
usual hZ and hA processes, while in CP-violating models, the H2Z and H1H3 signals are 
taken into account in addition to the dominant HiZ and H1H2 channels (the two other 
processes, H3Z and H2H3 being out of reach). 

3 Tools for the statistical analysis 

When scanning over the parameter space of a model, confidence levels are computed 
at each point to test the compatibility of data with the hypothesis of background only 
and with that of background plus signal as expected from the model. Throughout this 
section, the notations h, H and A must be understood as generic notations for the three 
neutral Higgs bosons of any type of MSSM scenario. 

3.1 Confidence level definitions and calculations 

The confidence levels are calculated using a modified frequentist technique based on 
the extended maximum hkelihood ratio [24] which has also been adopted by the LEP 
Higgs working group. The basis of the calculation is the hkelihood ratio test-statistic, Q: 

lnQ=-5 + ^ln^^ 

where S is the total signal expected and Sj and are the signal and background densities 
for event i. These densities are constructed using either expected rates only or also 
additional discriminant information, which can be one- or two-dimensional. Table 1 
presents the level of discriminant information for each channel: LEPl results rely on 
rates only, while LEP 2 results mix channels without or with discriminant information. 
As an example, in neutral Higgs boson channels with discriminant information, the first 
variable is the reconstructed Higgs boson mass in the hZ analyses and the sum of the 
reconstructed h and A masses in the hA analyses, while the second variable, if any, 
is channel-dependent, as specified in the references listed in the Table. Charged Higgs 
analyses use discriminant information in a similar way [23]. The searches for Higgs bosons 
decaying hadronically encompass analyses without or with Id discriminant information 
together with analyses whose selections vary with the mass hypothesis [22] . 
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The observed value of Q is compared with the expected Probabihty Density Func- 
tions (PDFs) for Q, which are built using Monte Carlo sampling under the assumptions 
that background processes only or that both signal and background are present. The 
confidence levels CLi, and CLg+b are their integrals from — oo to the observed value of 
Q. Systematic uncertainties in the rates of signal or background events are taken into 
account in the calculation of the PDFs for Q by randomly varying the expected rates 
while generating the distribution [25], which has the effect of broadening the expected Q 
distribution and therefore making extreme events seem more probable. 

CLb is the probability of obtaining a result as background-like or more so than the one 
observed if the background hypothesis is correct. Similarly, the confidence level for the 
hypothesis that both signal and background are present, CLg+b, is the probability, in this 
hypothesis, to obtain more background-hke results than those observed. The quantity CLg 
is defined as the ratio of these two probabilities, CLg+b/CLb. It is not a true confidence 
level, but a conservative pseudo-confidence level for the signal hypothesis. All exclusions 
discussed hereafter use CLg and require it to be 5% for an exclusion confidence of 95%. 
As using CLg instead of CLg+b is conservative, the rate of fake exclusions is ensured to 
be below 5% when CLg is equal to 5%. 

3.2 Estimation of expected signal and background densities 

The expected signal and background densities, which are required to check the con- 
sistency of the data with the background and signal processes have two components: the 
overall normalization which sets the expected rates and the Probability Density Functions 
(PDF) of the additional discriminant information, if any. 

The expected background and signal rates were calculated from the number of simu- 
lated events passing the cuts. For the signal the efficiencies derived from simulations at 
given mass points had to be interpolated to estimate efficiencies at Higgs boson masses 
which were not simulated. In most cases this was done using one polynomial or if neces- 
sary two polynomials, one to describe the slow rise, and a second to handle the kinematic 
cut-off, which can be much more abrupt. For the cases where two signal masses must be 
allowed, a two-dimensional parameterization was used. 

The shapes of the PDFs were derived using histograms which arc taken from the 
simulated events. In the case of two-dimensional PDFs these distributions were smoothed 
using a two-dimensional kernel, which consists of a Gaussian distribution with a small 
component of a longer tail [26]. The global covariance of the distribution was used to 
determine the relative scale factors of the two axes. The width of the kernel varied from 
point to point, such that the statistical error on the estimated background processes 
was constant at 20%. Finally multiphcative correction factors (each a one-dimensional 
distribution for one of the two dimensions of the PDF) were derived such that when 
projected onto either axis the PDF has the same distribution as would have been observed 
if it had been projected onto the axis first and then smoothed. This makes better use of 
the simulation statistics if there are features which are essentially one-dimensional, such 
as mass peaks. The error parameter fixed to 20% was an important choice. It was set by 
dividing the background simulation into two subsamples, generating a PDF with one and 
using the other to test for over-training by calculating the CLb obtained from simulation 
of background events. This should be 0.5 if the results are not to be biased, and a value 
of 20% for the error gave the closest approximation to 0.5 in all channels. Examples of 
smoothed two-dimensional PDFs are given in Fig. 3. 
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The signal simulations were made at fixed Higgs boson masses, but in order to test a 
continuous range of masses, interpolation software [27] was used to create signal PDFs at 
arbitrary masses. In the last year of operation, LEP energy was varied continuously while 
simulations were made at fixed beam energies. The same interpolation software was used 
to create signal and background PDFs at the correct ccntrc-of-mass energies [1]. The 
interpolation was done by linearly interpolating the cumulative distributions taking as a 
parameter the signal mass or the centre-of-mass energy. The procedure has been tested 
over ranges up to 40 GcV/c^ in mass while the actual shifts in the simulations were up to 
0.3 GeV in ^/s, and 5 GeV/c^ in mass for the hZ signals overall, but less than 0.5 GeV/c^ 
for Higgs boson masses between 113.5 and 116.5 GeV/c^. For the hA channels, the actual 
shifts were 5 GeV/c^ in either mass for Higgs boson masses between 80 and 95 GeV/c^ 
and up to 20 GeV/c^ elsewhere. Comparisons of simulated and interpolated distributions 
for a given mass were made in all channels and showed good agreement. 

3.3 The case of non-independent channels 

When combining the results in all channels to derive confidence levels, only inde- 
pendent channels must be included, which requires some special treatment for a few 
non-independent cases. 



a) different signals - one analysis with mass hypothesis-independent selections 



analysis 


(GeV) 


signals added 


ref. 


ffh four-b 


91 


(bbh ^ bbbb), (bbA bbbb), (hA ^ bbbb) 


[21] 


flh bbr+r" 


91 


(bbh ^ bbr+T"), (bbA ^ bbr+r"), (hA -> bbr+r") 


[21] 


ffh four-T 


91 


(r^r~h T^T~ T^T~), {t^t~ A t'^t~t'^t~) 


[21] 


wqq 


161-172 


(h -> qq) (Z -> i^u ), (h AA) {Z ^ i^u ) 


[17] 




189-208 


(h ^ r+r-) (Z ^ qq), (h ^ qq) (Z ^ r+r"), (hA ^ r+r-qq) 


[19,20,1] 


hZ four- jet 


161-183 


(h ^ qq) (Z ^ qq), (h ^ AA) (Z ^ qq) 


[17,18] 


hZ four- jet 


192-208 


(h ^ qq) (Z ^ qq), (h ^ AA) (Z ^ qq), (hA -> bbbb) 


[20,1] 


hA four-jet 


161-172 


(hA ^ bbbb), (h^ AA) A 


[17] 


hA four- jet 


192-208 


(hA ^ bbbb), (h ^ qq) (Z ^ qq) 


[20,1] 


four-b 


189-208 


(h -> AA)A, (h -> AA)Z, (hA -» bbbb) 


[21] 


b) different analyses - one final state 


final state 


^/i (GeV) 


competing analyses 


ref. 


four-jet 


91 


bbbb, bbcc 


[14] 


multi-jet 


91 


three and four- jet analyses 


[21] 




192-208 


low mass and high mass hZ analyses 


[1] 




189-208 


low mass and high mass flavour-blind analyses 


[22] 




189-208 


hZ and flavour-blind analyses 


[19,20,1,22] 


llqq,l=e,/i 


189-208 


hZ and flavour-blind analyses 


[19,20,1,22] 


four-jet 


192-208 


low mass and high mass hZ analyses 


[1] 




189-208 


low mass and high mass hZ fiavour-blind analyses 


[22] 




189-208 


three and four-jet hA flavour-blind analyses 


[22] 




189-208 


cscs and WAWA analyses 


[23] 




189-208 


hZ, hA, four-b, flavour-blind, cscs and W*AW*A analyses 


[19,20,1,21-23] 


TV jet jet 


189-208 


csTi>T and W*ATi>r analyses 


[23] 



Table 2: a) list of signals from the two lightest Higgs bosons h and A treated by a single 
analysis: the signal expectations are combined (rates added, PDFs summed with weights 
according to the rates) prior to the confidence level calculations, b) list of different 
analyses of the same final state: only one analysis is selected at each point in the scans, 
based on the best expected performance for exclusion. In this Table, h and A denote any 
neutral Higgs boson allowed to be produced in the indicated production processes. 
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The first case is that of different signals covered by the same analysis. The treatment 
of this depended upon whether the analyses were themselves independent of the mass 
hypothesis for the Higgs bosons. The set of search channels (see Table 1) contains mostly 
analyses of this kind. In that case, all signals selected by one analysis were combined 
into one global channel prior to the confidence level computation. Expected rates were 
added together and PDFs were summed with weights given by the expected rates of the 
individual signals. As an illustration, Table 2- a gives the hst of these signals and analyses 
on the example of the production of the lightest Higgs bosons, h and A, through the hZ 
and hA processes. When extending the combination to the third Higgs boson, H, the 
same procedure was followed, first for the various signals from that boson in the HZ and 
HA processes, and then to combine hZ and HZ signals or hA and HA signals. The PDF 
combination in such a case is illustrated in Fig. 3. 

A different procedure was applied in the case of different signals covered by the same 
analysis whose selections do depend on the mass hypothesis, as most searches of Ref. [22] 
do. Different signals are covered by these analyses only when including signals from the 
third Higgs boson, H. In that case, in each analysis only one signal (from either h or H) 
was selected at each point in the scanned parameter space and at each centre-of-mass 
energy, on the basis of the smallest expected CLg from experiments with no signal (that 
is, on the basis of the strongest average exclusion if no signal is present). 

The second case of non-independent channels is that of a large overlap in the events 
selected by different analyses sensitive to the same final state. The list of such analyses 
and final states is detailed in Table 2-b. Again, for each final state, only that analysis 
with the strongest expected exclusion power was retained at each test point. This is not 
optimal but ensures that the channels which are then combined in the global confidence 
level computations are independent. 

When the two cases just described (different signals covered by one analysis, different 
analyses sensitive to the same final state) were present simultaneously, the signal addition 
was performed before the final analysis selection. Then if that step involved more than 
two analyses, the final selection was made in successive iterations. To quote the four- 
jet final state as an example, at energies above 190 GeV, the total hZ and hA signals 
were first computed in each of the three four-jet analyses of Ref. [1] and in the four-b 
analysis of Ref. [21]. This summed three signals in the low and high mass hZ dedicated 
four-jet analyses ((h qq) (Z — > qq), (h — > AA) (Z — > qq) and hA — > bbbb), two 
signals in the hA dedicated four-jet analysis (hA — > bbbb and (h — > qq) (Z qq)) and 
three signals in the four-b analysis ((h AA)A, (h AA)(Z — >■ qq) and hA bbbb). 
The signals due the third Higgs boson, H, were computed in the same way and added 
to those from the h boson. Then, a choice was made between the low and high mass 
hZ dedicated four-jet analyses. The result of this selection was compared with the hA 
dedicated four-jet analysis, and the best of these was confronted with the four-b analysis. 
A choice was made between the remaining analysis and the best between the various 
fiavour-blind multi-jet analyses, that is the low mass and high mass hZ dedicated fiavour- 
bhnd analyses, and the three and four-jet hA dedicated flavour-bhnd analyses [22]. As 
multi-jet flavour-blind analyses use mass-hypothesis dependent criteria, selecting the best 
one implied also a choice between the h and H signals for each of them. The analysis 
retained was finally compared with the result of the selection between the two charged 
Higgs multi-jet analyses, the cscs and WAWA dedicated analyses [23]. 
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h|^|a, analysis at a/s=206.6 GeV 




Figure 3: An example of two-dimensional PDFs from the analysis of the hZ qq fi^fJ^^ 
channel at -^i = 206.6 GeV [1]. The first discriminant variable is built from the recon- 
structed Higgs boson mass while the second is the event b-tagging variable. Top, left: 
PDF for a hZ signal with rrih = 102 GeV/c^. Top, right: PDF for a HZ signal with 
mn = 115 GeV/c^. Bottom: PDF expected from the occurrence of both signals in a sce- 
nario where the expectations for the two signals are similar (cross-sections 32 and 42 fb, 
branching fractions into bb 92% and 91%, selection efficiencies 69% and 66% for hZ and 
HZ, respectively) leading to a double peak in the combined PDF. 
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4 The CP-conserving MSSM scenarios 

In most of the parameter space of the CP-conserving MSSM scenarios, only hZ and 
hA productions are kinematically possible at LEP energies. These processes have com- 
plementary cross-sections since the hZZ and hAZ couplings are proportional to sin(Q; — f3) 
and cos(a — P), respectively, where tan/3 is the ratio of the doublet vacuum expectation 
values and a is the Higgs doublet mixing angle which enters the definition of the two CP- 
even Higgs eigenstates as a mixture of the real, neutral components of the initial Higgs 
field doublets [2,28]. If kinematically allowed, hZ production dominates at low tan/5 or 
at large mA, while in the rest of the parameter space, it is suppressed with respect to 
hA pair-production. The third neutral Higgs boson, H, in some scenarios and in limited 
regions of the parameter space, is light enough and can be produced with a large HZ or 
HA cross-section. As the HZZ coupling is proportional to cos(a — /3), and the HAZ one 
is proportional to sin(a — f3), HZ production, when allowed by kinematics, plays a role 
at large tan/3, and HA production at low tan/5. Similarily, charged Higgs bosons kine- 
matically accessible at LEP2 energies are predicted in hmited regions of the parameter 
space, typically when A is light, whatever tan/3. The minimal value of the mass of such 
charged Higgs bosons is 60 GcV/c^ in the scenarios under study. The coverage of the 
region of the MSSM parameter space kinematically accessible at LEP is then assured 
primarily by the hZ and hA searches, with the help of the HZ, HA and to a lesser extent 
H+H~ channels. 

At tree level, the production cross-sections and the Higgs branching fractions in the 
MSSM depend on two free parameters, usually chosen as tan/5 and one Higgs boson 
mass, or, alternatively, two Higgs boson masses, e.g. uia and nih- Radiative corrections 
introduce additional parameters related to supersymmetry breaking [2,28]. Hereafter, the 
usual assumption that some of them are equal at a given energy scale is made: hence, 
the SU(2) and U(l) gaugino mass parameters are assumed to be unified at the so-called 
GUT scale, while the sfermion mass parameters or the squark trilinear couplings are 
taken to be equal at the EW scale. Within these assumptions, the parameters beyond 
tree level are: the top quark mass, the Higgs mixing parameter, /i, which defines the 
Higgsino mass parameter at the EW scale, the common sfermion mass parameter at the 
EW scale, Mgusy, the SU(2) gaugino mass parameter at the EW scale, M2, the gluino 
mass, rrig, and the common squark trilinear coupling at the EW scale, A. The U(l) 
gaugino mass term at the EW scale. Mi, is related to M2 through the GUT relation 
Ml — (5/3)tan^^vKM2. The radiative corrections affect the Higgs boson masses and 
couplings, with the largest contributions arising from loops involving the third generation 
quarks and squarks (top/stop and, at large values of tan/3, bottom/sbottom). As an 
example, the h boson mass, which is below that of the Z boson at tree level, increases 
by a few tens of GeV/c^ in some regions of the MSSM parameter space due to radiative 
corrections. 

4.1 The benchmark scenarios 

In the following, eight benchmark scenarios arc considered, as suggested in Ref. [29]. 
The values of their underlying parameters are quoted in Table 3. The first three scenarios 
are those usually studied at LEP. They have been proposed to test the sensitivity of 
LEP to Higgs bosons with either masses close to the kinematic limit or decays difficult 
to detect. Similarly, the five other scenarios are aimed at testing the sensitivity of the 
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scenario 


-^susy 

(GeV/c'') 


M2 

[GeV/c^) 


rUg ji Xt 
(GeV/c^) (GeV/c^) (GeV/c^) 




1000 


200 


800 


-200 


2 -^susy 


no mixing 


iOOO 


zOO 


800 


-200 





large 


400 


400 


200 


1000 


-300 


m^^, // > 


1000 


200 


800 


200 


2 -^susy 


m{f^^, /i > 0, < 


1000 


200 


800 


200 


-2 A^susy 


no mixing, > 0, large Mgusy 


2000 


200 


800 


200 





gluophobic 


350 


300 


500 


300 


-750 


small a 


800 


500 


500 


2.5 Mgusy 


-1100 



Table 3: Values of the underlying parameters for the eight representative MSSM scenarios 
scanned in this paper. Note that Xt — A — //cot/3. These scenarios have been studied 
for several values of the top quark mass, mtop = 169.2, 174.3, 179.4 and 183.0 GeV/c^. 



Higgs boson searches at hadron colliders. It is thus interesting to establish the LEP 
constraints in such models too. 

The first two scenarios, called the m™^ scenario and the no mixing scenario, differ 
only by the value of Xt = A — fi cot f3, the parameter which controls the mixing in the 
stop sector (through the product rritopXt). This parameter has the largest impact on 
the mass of the h boson. The m™^ scenario leads to the maximum possible h mass as 
a function of tan/3. The no mixing scenario is its counterpart with vanishing mixing, 
leading to theoretical upper bounds on nih which are at least 15 GeV/c^ lower than in 
the m™^^ scheme. 

The third scenario is called the large /j scenario to account for a large, positive value of 
11. As a consequence of the low value of Mgusy and the moderate mixing in the stop sector, 
this scenario predicts at least one CP-even Higgs boson with a mass within kinematic 
reach at LEP2 in each point of the MSSM parameter space. However, there are regions 
for which detecting such a Higgs boson is difficult because of vanishing branching fractions 
into b-quarks. The values chosen for /i and Xt are indeed such that, in these regions, 
radiative corrections lead to suppressed couplings to b-quarks for one or the other CP- 
even Higgs boson. The dominant decays in these regions being still into hadrons, the 
main analysis channels suffer from large backgrounds. This scenario was designed to 
test the sensitivity of LEP through analyses that could not benefit from the b-tagging 
capabilities of the experiments. 

Among the five other benchmark scenarios, three are variants of the m™^ and no 
mixing scenarios. The sign of fi and that of the mixing parameter have been reversed in 
the two scenarios derived from the LEP m™^^ scenario. The changes in the Higgs boson 
mass spectrum and properties are small. The sign of n has been reversed and the value 
of Mgusy has been doubled in the scenario derived from the no mixing scenario of LEP. 
The higher Msusy scale leads to a few GeV/c^ increase of the theoretical upper bound 
on TTT-h. The last two scenarios have been proposed to test potentially difficult cases 
for the searches at hadron colliders. Hence, the gluophobic scenario presents regions 
where the main production channel at the LHC, gluon fusion, is suppressed due to 
cancellations between the top quark and stop quark loops in the production process. 
Finally, in the small a scenario, important decay channels at the Tcvatron and at the 
LHC, h — > bb and h t~^t~, are suppressed at large tan/3 and moderate mA- In these 
regions, the radiatively corrected mixing angle a is low, resulting in suppressed couplings 
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scenario 


mtop (GeV/c^) 
169.2 174.3 179.4 183.0 


no mixing 
large 


128.2 132.9 138.6 142.7 

110Q "lice 1100 1 on o 
iiz.o iio.o llo.z iZU.O 

106.1 108.0 110.1 111.6 


m^^, > 

mjf*^^, ^ > 0, < 

no mixing, n > 0, large Mgusy 

gluophobic 

small a 


128.4 134.1 140.1 144.3 

124.5 128.8 134.3 138.2 
117.0 120.2 123.7 126.3 
115.7 118.8 122.0 124.4 
118.5 122.2 126.2 129.1 



Table 4: Maximal value of mh (in GeV/c^) in the eight benchmark MSSM scenarios 
studied in this paper, as a function of mtop- Radiative corrections include all dominant 
second-order loop terms [30]. The maximum value of mh corresponds approximately to 
the minimum value of the mass of the third Higgs boson, H. Bold values indicate scenarios 
where this boson is kinematically accessible at LEP. 

of the ligthest CP-even Higgs boson to down-type fermions since these couplings are 
proportional to — sin a/cos p. 

In all scenarios, the radiative corrections have been computed in the Feynman- 
diagrammatic approach with all dominant two-loop order terms included, using version 
2.0 of the FeynHiggs code [30]. As a first illustration of the different scenarios. Table 4 
gives the maximum value of mh allowed by theory in each of them, for the four values 
of mtop studied in this paper. At a given m^p value, the three m™*^^ scenarios give the 
highest upper bounds on mh, the positive ii scenario leading to the maximal value. The 
large /j, scenario presents the lowest upper bound, followed in increasing order by the no 
mixing scenario, the gluophobic one, the no mixing scenario with positive /i and the small 
a scheme. The maximum value of mh increases significantly with mtop- The effect is most 
important in the three m™^^ scenarios, and is much smaller in the others, especially in 
the large ^ scheme. It must be noted that the maximum value of mh corresponds approx- 
imately to the minimum value of mn in regions of large HZZ couplings (see Fig.4). Thus, 
there are a few scenarios where the H signal is expected to contribute to the experimental 
sensitivity. These are indicated in bold characters in Table 4. 

To illustrate further the Higgs boson phenomenology at LEP, a few properties are 
compared in Fig. 4 in the case of the no mixing and m™*^^ scenarios for a top quark mass 
of 174.3 GeV/c^. The figures showing masses and cross-sections underhne the importance 
of the signal from the heavy scalar, H, which can be kinematically accessible at LEP 2 
energies with a large HZ production cross-section at large tan/3 and moderate tua, up 
to about 100 GeV/c^. The width curves demonstrate that, at large tan/3, neutral Higgs 
bosons can have a width exceeding the experimental resolution which is of the order of 
1 to 3 GeV/c^ depending on the search channel. At moderate mA, this affects the h 
and A bosons and thus the hA production mode, but not the HZ one. At large mA, 
width effects become negligible for the h boson so that the hZ production mode, which 
is the only possible dominant mode in that region, is not affected. The figures showing 
branching fractions compare the no mixing and m™^ scenarios at low tan/3. In both 
scenarios, the h branching fraction into bb decreases to the profit of that into AA at very 
low tan/3, but the residual bb branching ratio is significantly higher in the m™'^'' scenario. 
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Figure 4: Properties of the three neutral Higgs bosons of the CP-conserving MSSM in 
the no mixing and m™*^^ scenarios with rritop = 174.3 GeV/c^. Top: H and h masses and 
H, h and A production cross-sections at ^/s = 206 GeV, at various tan P values. Middle: 
h, H and A widths as a function of ttia and tan jS. Bottom: h branching fractions as 
a function of mA at low to moderate values of tan/5. Decays into bb (solid lines) and 
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Finally, it should be noted that our previous MSSM interpretation of Ref. [1] relied 
on partial two-loop order radiative corrections [31]. In the present paper, these have been 
updated to include all dominant two- loop order corrections [30] . This leads to significant 
changes in the Higgs boson masses and properties. The main effect is an increase of the 
maximum (resp. minimum) allowed value of the h (resp. H) boson mass at fixed tan /?. 
As a consequence, the experimental sensitivity in tan/3 and that in mn are expected 
to decrease. A review of the changes induced by the more complete corrections on the 
experimental sensitivity of DELPHI is given in Ref. [32,33] in the framework of the three 
LEP scenarios, keeping identical experimental inputs. 

4.2 Scan procedure 

In each scenario, a scan was performed over the MSSM parameters tan/9 and mj^. 
The range in uia spans from 0.02 GeV/c^ up to 1 TeV/c^. Values of mA leading to 
unphysical negative mass squared values were removed from the scans. Such points are 
rather rare, except in the large n, gluophobic and small a scenarios (see section 5). The 
range in tan (3 extends from the minimal value allowed in each scenario^ up to 50, a value 
chosen in the vicinity of the ratio of the top- and b-quark masses, above which the Higgs- 
bottom Yukawa coupling is expected to become unreliable (see e.g. [2]). The scan steps 
were 1 GeV/c^ in ttia and 0.1 in tan/3 in the regions where mh varies rapidly with these 
parameters. At low ttia, where the decay modes change rapidly with the Higgs boson 
mass, values tested were 0.02, 0.1, 0.25, 0.5, 1.5 and 3 GeV/c^. 

At each point of the parameter space, the neutral and charged Higgs cross-sections and 
their branching fractions were taken from databases provided by the LEP Higgs working 
group [7] , on the basis of the theoretical calculations in Ref. [30] , completed by that in 
Ref. [34] for the charged Higgs boson branching fractions. The signals from the third Higgs 
boson, H, were included in the channel combination at each point where mn was found to 
be below 120 GeV/c^, the ultimate sensitivity of LEP. The signal expectations in each 
channel were then derived from the theoretical cross-sections and branching fractions, 
the experimental luminosity and the efficiencies. If necessary, a correction was applied 
to account for different branching fractions of the Higgs bosons between the test point 
and the simulation (e.g. for the hZ and HZ processes, the simulation was done in the SM 
framework) . 

As stated in the previous section, neutral Higgs bosons can have non-neghgible widths 
at large tan/3 when toa is above a few tens of GeV/c^. In this region, the experimental 

sensitivity is dominated by the LEP2 hA analyses dedicated to standard MSSM final 
states. To account for width effects in these channels, efficiencies derived from simula- 
tions with h and A widths below 1 GeV/c^ (see e.g. Ref. [1]) were applied for tan/3 < 
30 only. Above that value, efficiencies were linearly interpolated in tan/3 between the 
efficiencies from these simulations and those from simulations at tan/3 = 50 where the 
Higgs boson widths exceed the experimental resolution. As the Higgs boson widths grow 
approximately linearly with tan (3 above 30, a linear interpolation is valid. The same holds 
for the discriminant information, for which the same interpolation software was used as 
discussed in section 3.2 for the PDF interpolation in mass or centre-of-mass energy. The 
effect of the Higgs boson widths on the PDFs of the hA signals and the interpolation in 
tan /3 of these PDFs arc illustrated in Fig. 5. Note that the hZ and HZ channels at large 
tan j3 are not affected by such an effect since in most of the regions where they possibly 
contribute, their widths are below the experimental resolution, as shown in Fig. 4. 

"^Thc minimal value of tan /3 is around 0.7 in the large ^ scenario and in the no mixing scenario with positive /i and 0.4 
in all other schemes. Lower tan /3 values give rise to unphysical negative mass squared values in the Higgs sector. 
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Figure 5: Two-dimensional PDFs used in the analysis of the hA bbbb channel at 
\fs = 206.6 GeV [1]. The first discriminant variable is the sum of the reconstructed 
Higgs boson masses while the second is a neural network output. Top, left: PDF for a 
hA signal with rriA = '^h = 90 GeV/c^ and h and A widths below 1 GeV/c^. Top, right: 
PDF for a hA signal with ttia = ?7^h = 90 GeV/c^ and tan/3 = 50. The Higgs boson 
widths in that case are 5 and 9 GeV/c^ for A and h, respectively. Bottom: PDF linearily 
interpolated in tan at a value of 37. 
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5 Results in CP-conserving MSSM scenarios 

The regions of the MSSM parameter space excluded at 95% CL or more by combining 

the searches of Table 1 are hereafter discussed in turn for each scenario. The exclusion is 
dominated by the searches for neutral Higgs bosons in standard MSSM final states. The 
searches for neutral Higgs bosons decaying into hadrons of any flavour and the charged 
Higgs boson searches complete the exclusion in restricted regions of the parameter space. 
In addition, the limit on the Z partial width that would be due to new physics [35], 
-pnew ^ Q Q MgY/c^ is used as an external constraint on the hA process at LEPl. A 
detailed account of the impact of these auxiliary constraints can be found in Ref. [33,36]. 

5.1 The m™^^ scenario 

The excluded regions in the m™*^^ scenario are presented in the (mt, tan/?), (mA, 
tan P) and (mh, planes in Fig. 6 for a top mass value of 174.3 GeV/c^. Basically, the 
exclusion is made by the results in the hZ (hA) channels in the low (large) tan /3 region 
while they both contribute at intermediate values. The searches for the heavy scalar, H, 
brings no additional sensitivity since H is not kinematically accessible in this scenario 
(see Table 4). The above results establish the following 95% CL lower limits on mh and 
niA for mtop = 174.3 GeV/c^: 

mh > 89.7 GeV/c^ tua > 90.4 GeV/c^ 

for any value of tan (3 between 0.4 and 50. The expected median limits are 90.6 GeV/c^ for 
mh and 90.8 GeV/c^ for mA- The observed limit in mA (mh) is reached at tan P around 20 
(10), in a region where both the hZ and hA processes contribute. For mtop = 174.3 GeV/c^ 
the range in tan/? between 0.7 and 1.9 (expected [0.7-1.9]) is excluded for any value of 
mA between 0.02 and 1000 GeV/c^. These limits and exclusions, as well as those for all 
the CP-conserving scenarios, are summarized in Table 5. 

The mtop dependence of the above limits was also studied, as reported in Table 5. The 
mass limits remain unchanged when varying mtop, for mh is insensitive to mtop in the 
region of large tan/? and intermediate mA where the limits are set. On the other hand, 
the excluded range in tan/3 is governed by the maximal value of mh, which is reached 
at large mA where mh is very sensitive to mtop, as illustrated in the top left-hand plot 
in Fig. 6: hence the variation of the limits in tan/? as reported in Table 5 and Fig. 14. 
An exclusion in tan/? exists for a top mass up to 179.4 GeV/c^ which is about three 
standard deviations higher than the current average mtop measurement. The exclusion 
would vanish for a top mass as high as 183.0 GeV/c^. 

5.2 The scenario but with /i positive and either sign for Xt 

The excluded regions for a top mass value of 174.3 GeV/c^ are presented in Fig. 7 for 
the m™^ scenario with positive /i, keeping Xt positive as in the original m^^^ scenario, 
and in Fig. 8 for the m™^ scenario with positive and negative Xf. The results are 
quite similar to those in the original m™'^^ scenario. Mass limits arc within 200 MeV/c^ 
of those in the previous section and do not vary significantly with mtop, as reported in 
Table 5. 

To compare observed and median limits, the 95% CL lower limits on mh and mA in 
the m™^ scenario with positive fi for mtop — 174.3 GeV/c^ are: 
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Figure 6: MSSM rn^^^ scenario for a top mass of 174.3 GeV/c^: regions excluded at 
95% CL by combining the results of the Higgs boson searches in the whole DELPHI 
data sample (light-grey). The dashed curves show the median expected limits. The 
medium-grey areas are the regions not allowed by theory. The dash-dotted lines in the 
top left-hand plot are the theoretical upper bounds for a top mass of 169.2, 179.4 and 
183.0 GeV/c2 (from left to right). 
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Figure 7: MSSM rri^^^ scenario with positive /i for a top mass of 174.3 GeV/c^: regions 
excluded at 95% CL by combining the results of the Higgs boson searches in the whole 
DELPHI data sample (light-grey). The dashed curves show the median expected limits. 
The medium-grey areas are the regions not allowed by theory. The dash-dotted lines in 
the top left-hand plot are the theoretical upper bounds for a top mass of 169.2, 179.4 and 
183.0 GeV/c2 (from left to right). 
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Figure 8: MSSM m™^^ scenario with positive /x and negative Xt for a top mass of 
174.3 GeV/c^: regions excluded at 95% CL by combining the results of the Higgs boson 
searches in the whole DELPHI data sample (light-grey). The dashed curves show the 
median expected limits. The medium-grey areas are the regions not allowed by theory. 
The dash-dotted lines in the top left-hand plot are the theoretical upper bounds for a top 
mass of 169.2, 179.4 and 183.0 GeV/c^ (from left to right). 
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mh > 89.6 GeV/c^ 



rriA > 90.3 GeV/c^ 



for any value of tan/3 between 0.4 and 50. The expected median limits are 90.3 GeV/c^ 
for rrih and 90.4 GeV/c^ for ijia. The 95% CL lower limits on mh and ijia in the m™^ 
scenario with positive /i and negative Xt for mtop = 174.3 GeV/c^ are: 



for any value of tan/3 between 0.4 and 50. The expected median limits are 90.4 GeV/c^ 
for rrih and 90.6 GeV/c^ for mA. 

The excluded ranges in tan /3 are different in the three m™^ scenarios, since they have 
different theoretical upper bounds on mt. For mtop = 174.3 GeV/c^ the excluded range 
in the m™*^^ scenario with positive lies between 0.7 and 2.0 (expected [0.8-2.0]), while 
in the m^^^ scenario with positive and negative Xt it spans from 0.6 to 2.5 (expected 
[0.6-2.4]). These hmits are valid for any value of mA between 0.02 and 1000 GeV/c^. 
Note that despite the higher maximal value of mh in the m^^'^ scenario with positive /i, 
the most conservative limits in tan f3 are still derived in the original m]f scenario (see 
Section 5.1), reflecting the differences in the theoretical upper bounds at tan/3 around 1 
(see top left-hand plots in Fig. 6, 7 and 8). The mtop dependence of the above limits is 
presented in Table 5 and Fig. 14. For a top mass as high as 183 GeV/c^, there would be 
no longer any exclusion in tan /3 in the m'^'^^ scenario with positive n, while there is still 
one in the scenario with positive /i and negative Xt due to the lower maximal value of 
mh in that scenario. 

5.3 The no mixing scenario 

The excluded regions in the no mixing scenario are presented in Fig. 9 for a top mass 
value of 174.3 GeV/c^. In this scenario, if the top is not too heavy, the heavy scalar, H, 
is kinematically accessible at large tan/3 and moderate mA, the region where the mass 
limits in mA and mh are set. Thus, allowing for its production increases the sensitivity 

of the searches. 

The zoom at low mA in the (mh, mA) projection shows that the direct searches leave 
three unexcluded regions below 12 GeV/c^ in mA. The thin strip along the theoretical 
lower bound on mh at very low mA (hardly visible in the figure) is excluded by the limit 
on the Z partial width that would be due to new physics [35], F"^*^" < 6.6 MeV/c^, which, 
when applied to the hA process, translates into an excluded region that encompasses 
that area. This is not the case for the two other unexcluded regions. These have tan/3 
below 1.0 and mh between 59 and 82 GeV/c^. In that region, mA is below the kinematic 
threshold mh = 2mA, the decay h — > A A opens and supplants the h — > bb mode, as can 
be seen in Fig. 4. Our LEP2 h — * AA searches, covering A masses above the cc threshold 
(see Table 1), have no sensitivity below 4 GeV/c^ in mA- Similarly, charged Higgs bosons, 
although kinematically accessible with a mass between 57 and 82 GeV/c^, have a large 
branching fraction into W*A in this region. As our charged Higgs boson searches in these 
channels assume mA above 12 GeV/c^ (see Table 1), the overall experimental sensitivity 
in these regions remains weak and no exclusion at 95% CL can be derived, in agreement 
with the expected performance. The largest value of CLg is 7% in the unexcluded region 
around 12 GeV/c^ in mA and 33% in the unexcluded hole below 4 GeV/c^. Note that the 
nearby region with mh from 82 GeV/c^ to the theoretical upper bound on mh is excluded 
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Figure 9: MSSM no mixing scenario for a top mass of 174.3 GeV/c^: regions excluded 
at 95% CL by combining the results of the Higgs boson searches in the whole DELPHI 
data sample (light-grey). Among the three unexcluded regions at low mA, the strip at 
low mh is fully excluded by the limit on the Z partial width that would be due to new 
physics [35]. The dashed curves show the median expected limits. The medium-grey 
areas are the regions not allowed by theory. The dash-dotted lines in the top left-hand 
plot are the theoretical upper bounds for a top mass of 169.2, 179.4 and 183.0 GeV/c^ 
(from left to right). 
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at 95% CL by the charged Higgs boson searches through their fermionic decays which 
dominate the W*A mode there. 

The above results estabhsh the following 95% CL lower hmits on rrih and uia for 
mtop = 174.3 GeV/c^: 

mh > 90.7 GeV/c2 tjia > 91.2 GeV/c^ 

for any value of tan/? between 1.0 and 50. The expected median limits are 91.1 GeV/c^ 
for both mh and rriA- The observed limits in ttia and mh are reached at tan/3 around 15, 
in a region where both the hZ and hA processes contribute. For mtop = 174.3 GeV/c^, 
two ranges in tan P are excluded for any value of uia between 0.02 and 1000 GeV/c^, the 
largest interval being between 1.0 and 9.7 (expected [0.9-7.7]). 

The iTT-top dependence of the above limits was studied, as shown in Table 5 and Fig. 14. 
In this scenario, both the mass limits and the excluded range in tan f3 change when varying 
mtop- Indeed, as already mentioned, the mass limits in ttia and mh rely on the searches for 
H, whose mass is very sensitive to m^p in the region where the limits are set. Similarly, 
the maximal value of mh, which governs the limits in tan/5, is reached at large rriA where 
mh is very sensitive to mtop (see Table 4). Note that for a top mass of 169 GeV/c^, 
mn decreases by 3 GeV/c^ in the region where the mass limits are set, making the H 
signal more within the sensitivity of LEP2: the whole parameter space of the no mixing 
scenario is then accessible and found to be excluded at 95% CL, apart from two holes 
at tan/? below 1.0, one at rriA around 12 GeV/c^, which is excluded at 92% CL, and a 
larger one below 4 GeV/c^, which is disfavoured at 69% CL only. 

5.4 The no mixing scenario but with positive /i and large Mgusy 

The excluded regions in the no mixing scenario with positive ^ and large Mgusy are 
presented in Fig. 10 for a top mass value of 174.3 GeV/c^. The larger Mgusy makes 
the impact of the H signal, and hence the exclusion limits, weaker than in the previous 
scenario. On the other hand, the results in the low mass region, at uia below 12 GeV/c^, 
are similar to those in the no mixing scenario. The direct searches leave a tiny unexcluded 
strip at low mh and very low itla which is excluded by the limit on T'"^". Three other 
regions, at mh between 56 and 72 GeV/c^, remain unexcluded even when charged Higgs 
boson searches are included, due to the large branching fraction into W*A decays, which 
are not covered by these searches at such low A masses. The holes around 8 and 12 GeV/c^ 
in TTiA are however excluded at 93% and 91% CL, respectively, while the larger area below 
4 GeV/c^ in ruA is disfavoured at 60% CL only. 

The above results establish the following 95% CL lower hmits on mh and niA for 
mtop = 174.3 GeV/c^: 

mh > 89.8 GeV/c^ ttia > 90.6 GeV/c^ 

for any value of tan/5 between 1.0 and 50. The expected median limits are 90.5 GeV/c^ 
for mh and 90.6 GeV/c^ itla- For mtop = 174.3 GeV/c^ the range in tan/3 between 1.0 and 

4.5 (expected [1.0-4.3]) is excluded for any value of uia between 0.02 and 1000 GeV/c^. 
The mtop dependence of the above limits is presented in Table 5 and Fig. 14. The 

mass limits vary only slightly with mtop, since in the region where these are set, mh is 
insensitive to mtop while mn, although sensitive to mtop, is very close to the kinematic 
limit. Contrary to the case of the no mixing scenario, the parameter space of this scenario 
does not become fully accessible for a top mass of 169 GeV/c^, due to too high an upper 
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Figure 10: MSSM no mixing scenario with positive /i and large Mgusy for a top mass of 

174.3 GeV/c^: regions excluded at 95% CL by combining the results of the Higgs boson 
searches in the whole DELPHI data sample (light-grey). Among the four unexcluded 
regions at low m^, the strip at low mh is fully excluded by the limit on the Z partial width 
that would be due to new physics [35]. The dashed curves show the median expected 
limits. The medium-grey areas are the regions not allowed by theory. The dash-dotted 
lines in the top left-hand plot are the theoretical upper bounds for a top mass of 169.2, 

179.4 and 183.0 GeV/c^ (from left to right). 
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(resp. lower) bound on rrih (resp. mn). The exclusion is thus much weaker than in the 
no mixing scheme but stronger than in the m™^ scenarios. 

5.5 The large fi scenario 

The excluded regions in the large fi scenario are presented in the (mh, tan/?) and (ttia, 
tan/3) planes in Fig. 11 for values of the top quark mass of 174.3 and 179.4 GeV/c^. In 
these figures, the contribution of the H signal and that of the searches for neutral Higgs 
bosons decaying into hadrons of any flavour are highlighted. 

A large fraction of the allowed domain is excluded by the searches for the h, A and H 
Higgs bosons into standard MSSM final states. In particular, since the theoretical upper 
bound on the h boson mass in this scenario is low (around 110.0 GeV/c^, see Table 4), 
the sensitivity of the hZ channels is high even at large tan /3, which explains why the 
excluded region reaches the theoretically forbidden area for large values of tan f3. As the 
value of the upper bound on mh is also the theoretical lower bound on mn at large tan f3, 
allowing for the production of H translates into a significant gain in exclusion, namely at 
tan P above 8. The searches for neutral Higgs bosons decaying into hadrons of any flavour 
bring an additional exclusion in regions left uncxcluded by the standard searches at tan /3 
around 14. At moderate mA, hZ and hA productions are low due to weak hZZ couplings 
for hZ and to kinematics for hA. On the other hand, HZ production is large but H is 
decoupled from bb. At larger uia, hA and HZ productions are kinematically forbidden, 
hZ production is large but the h— i^bb branching fraction vanishes. In both cases, the 
Higgs boson whose production is allowed (H or h) has a large branching fraction into 
hadrons and a mass close to the sensitivity of our searches for a neutral Higgs boson 
decaying into hadrons and fully coupled to the Z. This explains why these searches lead 
to an additional but only partial exclusion in these regions. Note that increasing the 
top quark mass from 174.3 to 179.4 GeV/c^ leads to a larger unexcludcd area. There 
are indeed more points with vanishing h or H branching fractions into bb and, as 
and mn increase with mtop, the impact of the searches for hadronically decaying Higgs 
bosons also becomes weaker. However, when combining the four LEP experiments, the 
sensitivity of these searches increases and becomes high enough to cover almost entirely 
these regions of vanishing branching fractions into bb [7]. 

Below 3 GeV/c^ in tua, the direct searches leave two unexcluded holes at tan j3 around 
1. The one at tan P above 1 is fully excluded by the limit on F"^^ for either value of mtop- 
The hole at tan /3 below 1 remains unexcluded. The largest value of CLg in this area is 
12% for mtop = 174.3 GeV/c and 6% for mtop = 179.4 GeV/c. 

The above results establish the following 95% CL lower limits on mh and mA for 
mtop = 174.3 GeV/c^: 

mh > 94.2 GeV/c^ mA > 96.6 GeV/c^ 

for any value of tan/? between 0.9 and 50. The expected median hmits are 90.3 GeV/c^ 
for mh and 92.8 GeV/c^ for mA. The observed limits in mA and mh are reached at 
tan (3 around 14, in a region where the hZ, HZ and hA processes contribute. For 
^top = 174.3 GeV/c^, two ranges in tan/? are excluded for any value of mA between 
0.02 and 1000 GeV/c^, the largest interval being between 0.9 and 13.7 (expected [0.9- 
12.9]). 

The mtop dependence of the above limits is presented in Table 5 and Fig. 14. Ex- 
cept for mtop = 174.3 GeV/c^, the mass limits vary only slightly with mtop and are in 
agreement with the expected ones. The difference at mtop = 174.3 GeV/c^ has been 



25 




20 40 60 80 100 120 140 20 40 60 80 100 120 140 

(GeV/c^) nij, (GeV/c^) 




50 100 150 200 250 300 350 400 450 500 50 100 150 200 250 300 350 400 450 500 

(GeV/c^) (GeV/c^) 



Figure 11: MSSM large fi scenario: regions excluded at 95% CL by combining the 
results of the Higgs boson searches in the whole DELPHI data sample (light-grey area 
and embedded domains in medium- and dark-grey). Results are shown for two values 
of the top mass, 174.3 and 179.4 GeV/c^. The domains embedded in the light-grey 
area at large tan/3 are excluded by the searches for the heavy scalar Higgs boson, H 
(medium-grey or green) and by the flavour-blind searches (dark-grey or dark-blue). Of 
the two unexcluded holes at low ttla, the one at tan/3 above 1 is excluded by the limit on 
the Z partial width [35] that would be due to new physics. The dashed curves show the 
median expected limits. The medium-grey areas with bold contours are the regions not 
allowed by theory. Note in particular the large region forbidden at low mj^ in the (mA, 
tan (3) projections, which is due to points leading to unphysical h masses. 
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traced back to the deficit in data with respect to background expectations which was 
observed in the fiavour-bhnd searches apphed to the Higgsstrahlung process [22] when 
testing masses above 100 GeV/c^, which corresponds to the range of mn values in the 
region where the mass limits are obtained in the large /i scenario. In this region, the set of 
independent channels which are selected to be statistically combined (sec Sec. 3.3) varies 
strongly from one top mass value to the other, due to still large H branching fractions 
into bb at mtop — 169.2 GeV/c^ and to mn values increasing with mtop (see Table 4). 
At mtop = 174.3 GeV/c^, the weight of the flavour-blind HZ searches is maximal and 
the deficit in data of these searches translates into a difference between the observed and 
median limits. 

5.6 The gluophobic scenario 

For the gluophobic scenario the excluded regions in the (mh, tan/5), (mA, tan/5) 
and {mil, ^a) planes arc presented in Fig. 12 for a top mass value of 174.3 GeV/c^. 
Although this scenario was designed to test Higgs boson searches at hadron colliders, 
with a phenomenology very different from that of LEP, results are similar to those 
derived in the previous scenarios. The exclusion is defined by the results in the hZ (hA) 
channels in the low (large) tan f3 region while they both contribute at intermediate values. 
The direct searches leave several unexcluded holes below 4 GeV/c^ in niA and at tan/5 
below 2, which are all excluded by the limit on F"^^. 

The above results establish the following 95% CL lower limits on mh and ttia for 
mtop = 174.3 GeV/c^: 

mh > 87.0 GeV/c^ mA > 92.9 GeV/c^ 

for any value of tan/3 between 0.4 and 50. The expected median limits are 87.0 GeV/c^ 
for mh and 93.0 GeV/c^ for mA- The observed limits in mA and mh are reached at tan/5 
around 50, in a region where only the hA process contributes. Contrary to the other 
scenarios, the h and A bosons are not degenerate in mass at large tan /5, which reflects 
in the significant difference between the h and A mass limits. For mtop = 174.3 GeV/c^, 
the range in tan/5 between 0.4 and 5.2 (expected [0.4-4.8]) is excluded for any value of 
mA between 0.02 and 1000 GeV/c^. 

The mtop dependence of the above limits is shown in Table 5 and Fig. 14. As already 
mentioned, the h and A bosons are not degenerate at large tan/5 and moderate mA, the 
region where the mass limits are set. As a consequence, the value of mh at fixed mA and 
tan /5 is observed to vary significantly with mtop in that region. This is the main reason 
for the variations of the mass limits with mtop, an additional effect being the variations 
of mn which is kinematically accessible at low mtop in this scenario (see Table 4). On the 
other hand, the variation of the excluded range in tan /5 is due, as in the other scenarios, 
to the change in the maximal value of mh which is very sensitive to mtop- 

5.7 The small a scenario 

The excluded regions in the small a scenario are presented in Fig. 13 for a top mass 
value of 174.3 GeV/c^. The small a scheme is the second example of a scenario aiming 
at testing potentially difficult cases for the Higgs boson searches at hadron colliders. As 
mentioned in section 4, this scenario presents regions of the parameter space where the h 
— > bb and h r'^r^ decays vanish, which could be a problem at LEP too. The results in 
Fig. 13, similar to those derived in the previous scenarios, show that this is not the case. 
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Figure 12: MSSM gluophobic scenario for a top mass of 174.3 GeV/c^: regions excluded 
at 95% CL by combining the results of the Higgs boson searches in the whole DELPHI 
data sample (light-grey). The unexcluded holes at low rriA are fully excluded by the 
limit on the Z partial width [35] that would be due to new physics (dark-grey). The 
dashed curves show the median expected limits. The medium-grey areas are the regions 
not allowed by theory. Note in particular the large forbidden region in the (mA, tan/3) 
projection, which is due to points leading to unphysical h masses. The dash-dotted lines 
in the top left-hand plot are the theoretical upper bounds for a top mass of 169.2, 179.4 
and 183.0 GeV/c^ (from left to right). 
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Figure 13: MSSM small a scenario for a top mass of 174.3 GeV/c^: regions excluded at 
95% CL by combining the results of the Higgs boson searches in the whole DELPHI data 
sample (light-grey). There is one unexcluded hole at low mh and tan/3 around 20 which 
is excluded by the limit on the Z partial width [35] that would be due to new physics 
(dark-grey). The dashed curves show the median expected limits. The medium-grey 
areas are the regions not allowed by theory. Note in particular the large forbidden region 
in the (mA, tan /3) projection, which is due to points leading to unphysical h masses. The 
dash-dotted lines in the top left-hand plot are the theoretical upper bounds for a top 
mass of 169.2, 179.4 and 183.0 GeV/c^ (from left to right). 
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At large tan/3, in the region accessible at LEP, the h — > bb branching fraction, although 
reduced, remains high enough (e.g. above 70% in the region where the mass limits are 
set) to ensure a good sensitivity. At low mh, the direct searches leave one unexcluded 
island that is fully excluded by the limit on T"*^™. 

The above results establish the following 95% CL lower limits on mh and mx for 
mtop = 174.3 GeV/c^: 

mh > 83.5 GeV/c^ mx > 95.8 GeV/c^ 

for any value of tan/3 between 0.4 and 50. The expected median hmits are 82.6 GeV/c^ 
for mh and 95.0 GeV/c^ for mA- The observed limits in mx and mh are reached at 

tan/3 around 50, in a region where only the hA process contributes. As in the previous 
scenario, the h and A bosons are not degenerate in mass at large tan /3, which reflects in 
the significant difference between the h and A mass limits. For m^p = 174.3 GeV/c^, 
the range in tan/3 between 0.4 and 4.0 (expected [0.5-3.9]) is excluded for any value of 
mA between 0.02 and 1000 GeV/c^. 
The 

^top dependence of the above limits is shown in Table 5 and Fig. 14. As in the 
previous scenario, the value of mh at fixed mA and tan /3 varies significantly with mtop in 
the region where the mass limits are set, which explains the variations of the latter. The 
H signal, being kinematically inaccessible for most values of mtop (see Table 4) plays no 
role in this scenario. Finally, the variation of the excluded range in tan /3 is due to the 
change in the maximal value of mh which is very sensitive to mtop- 

5.8 Summary 

The lower bounds in mass and excluded ranges in tan/3 obtained in the eight CP- 
conserving benchmark scenarios presented in the previous sections are summarized in 
Table 5. The variation with mtop of the excluded ranges in tan /3 is further illustrated in 
Fig. 14. All lower bounds in mass are at the 95% CL, as well as each individual (either 
lower or upper) bound in tan /3. In all scenarios, the radiative corrections on the Higgs 
boson masses and couplings have been computed in the Feynman-diagrammatic approach 
with all dominant two-loop order terms included, using version 2.0 of the FeynHiggs 
code [30]. 
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mop (GeV/c^) 




scenario 


limits 


169.2 


174.3 


179.4 


183.0 






89.7 


89.7 


89.7 


89.6 




niA 


90.4 


90.4 


90.4 


90.4 




tan P 


0.59 - 2.46 


72 - 1 96 


93 - 1 46 


none 


max 
'"'h 




89.6 


89.6 


89.5 


89.6 


/X > 


777, A 


90.3 


90.3 


90.3 


90.3 




tan /3 


0.59 - 2.61 


0.71 - 2.00 


0.87 - 1.54 


none 






89.6 


89.6 


89.5 


89.6 


u > X+ < 


777, A 


90.5 


90.4 


90.4 


90.4 




tan B 


53 - 3 20 


0.63 - 2.46 


0.72 - 1.96 


84 - 1 63 


no mixing 


777 V, 


112.8 


90.7 


90.0 


89.9 




777a 


1000. 


91.2 


90.8 


90.5 




tan B 


40 - 50 


0.40 - 9.70 


0.40 - 5.40 


0.40 - 4.40 




tan /9, 7r7A>12.0 


0.46 - 0.96 


46 - 96 


0.47 - 0.97 


0.47 - 0.97 


no mixing 


777,1, 


89.9 


89.8 


89.7 


89.8 


/X > 


777a 
' ' "A 


90.8 


90.6 


90.4 


90.3 


large Msusy 


tan P 


0.70 - 6.95 


0.70 - 4.55 


0.70 - 3.43 


0.70 - 2.97 




tan/3, 777a>12.0 


0.70 - 1.01 


0.70 - 1.01 


0.70 - 1.02 


0.70 - 1.01 


Large 


777h 


90.2 


94.2 


89.7 


89.3 




777A 


92.5 


96.6 


92.6 


92.5 




tan/3 


0.72 - 14.79 


0.72 - 13.68 


0.72 - 10.91 


0.72 - 10.63 




tan/3, 777a>2.4 


0.72 - 0.79 


0.75 - 0.85 


0.86 - 0.90 


none 


Gluophobic 


777h 


87.8 


87.0 


86.4 


86.2 




777A 


93.0 


92.9 


93.2 


93.5 




tan/3 


0.40 - 9.70 


0.42 - 5.22 


0.48 - 3.76 


0.51 - 3.19 


Small a 


777h 


84.3 


83.5 


82.5 


82.0 




777A 


95.0 


95.8 


96.5 


97.2 




tan /3 


0.40 - 5.97 


0.43 - 4.03 


0.52 - 3.12 


0.55 - 2.69 



Table 5: 95% CL lower bounds on 777h and 777a in GeV/c^ and excluded ranges in tan/3 
obtained in the different MSSM CP-conserving benchmark scenarios, as a function of 
mtop- Except for the two no mixing and the large fi scenarios, the exclusions in mass are 
valid for all values of tan/3 between 0.4 and 50, and the exclusions in tan/3 hold for all 
values of 777a between 0.02 and 1000 GeV/c^. In the three other scenarios, part of the 
interval in tan/3 is excluded only for rriA above a few GeV/c^ threshold: this sub-interval 
is indicated in a fourth line together with the threshold in ttla- As a consequence, the 
mass bounds in these scenarios are valid only for values of tan /3 outside the quoted 
sub- interval. 
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Figure 14: Variation with mtop of the ranges in tan/5 excluded by DELPHI in the 
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IJ, scenario. 
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6 The CP-violating MSSM scenarios 

In most of the parameter space of the CP- violating MSSM scenarios studied in this 
paper, only the two lightest neutral Higgs bosons, Hi and H2 are kinematically accessible 
at LEP energies. If their couplings to the Z boson are not strongly suppressed by CP- 
violation, the main production processes are the HiZ, H2Z and H1H2 processes, with HiZ 
dominating at low tan /3, H1H2 at large tan /5 and H2Z contributing over the whole range 
of tan (3 values allowed in each scenario. In restricted areas of the parameter space, the 
second pair-production process, H1H3, can add a non-negligible signal and has also been 
considered in the searches. On the other hand, in most scenarios, charged Higgs bosons 
have a mass above 100 GeV/c^, and thus have not been included. 

As already mentioned, CP violation in the MSSM Higgs sector is introduced through 
radiative corrections. Besides the two parameters used to define the scenarios at tree level, 
chosen as tan j3 and mH± , radiative corrections introduce additional parameters. As in the 
CP-conserving case, these are primarily mtop and the set of parameters related to super- 
symmetry breaking: fj,, Mg^sy, M^, rrig and A, as defined in section 4 [2,28]. In addition, 
CP violation introduces phases. The unification assumptions made for the supersym- 
mctry breaking parameters, and the global symmetries that govern the dimension-four 
operators of the MSSM Lagrangian, can be used to reduce the number of CP-violating 
phases to only two [37]. In the scenarios studied hereafter, these phases are taken as 
the phase of the gluino mass, arg(mg) and the phase of the common stop and sbottom 
trilinear coupling, arg(A). 



6.1 The benchmark scenarios 

The dominant CP-violating effects on the neutral Higgs boson masses and couplings 
to gauge bosons are proportional to 

mj^ Im(iiA) 

^ susy 

where v'^ is the quadratic sum of the vacuum expectation values of the two Higgs field 
doublets [38]. Sizeable effects are thus expected for moderate values of Mg^sy, large values 
of jji and phases arg(A) around 90°. A strong dependence on the value of mtop is also to 
be expected. 

Along these lines, Ref. [38] proposed a benchmark scenario with maximal CP- violation, 
the CPX scenario, as an appropriate scheme for direct searches at LEP and other col- 
liders. The values of its underlying parameters are quoted in Table 6. As expected from 
the above discussion, the value of M^^^y, a few hundred GeV/c^, is moderate, /i and 
\A\ take large values, 2 and 1 TeV/c^ respectively, and the CP- violating phase arg(yl) 
is set at 90°. Although the gluino- mass phase has a small impact on the CP-violating 
effects, these appear to be reinforced at 90° [37], a value which was thus retained for 
arg(mg). The values listed in Table 6 fulfill the existing constraints from measurements 
of the electron and neutron electric dipole moments, by making the first two generations 
of squarks sufficiently heavy, with masses above 1 TeV/c^. In the following, the CPX 
scenario has been studied for four values of the top quark mass, mtop — 169.2, 174.3, 
179.4 and 183.0 GeV/c^. 

In addition to the CPX scenario, a few variants have also been considered in order to 
study the dependence of the CP- violation effects on the values of phases, ^ and Msusy 
The values tested are quoted in Table 6. The two CP-violating phases, still taken to be 
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scenario 


M,u.y Ma 
(GeV/c^) (GeV/c^) 


\mg\ 
(GeV/c^) 


(GeV/c^) 


1^1 
(GeV/c^) 


arg(mg)=arg(A) 
(degrees) 


CPX 


500 


200 


1000 


2000 


1000 


90 


phase study 


500 


200 


1000 


2000 


1000 


0,30,60,135,180 


fj, study 


500 


200 


1000 


500,1000,4000 


1000 


90 


M,u,y =1 TeV/c^ 
Msusy =1 TeV/c^, scaled 


1000 


200 


1000 


2000 


1000 


90 


1000 


200 


2000 


4000 


2000 


90 



Table 6: Values of the underlying parameters for the representative CP-violating MSSM 
scenarios scanned in this paper, namely the CPX scenario and its ten variants. 

equal, were varied from to 180°, keeping all other parameters as in the CPX scenario. 
Values of fi below and above 2 TeV/c^ were studied in the same way. Finally, the value 
of Msusy was increased from 500 GcV/c^ to 1 TeV/c^, keeping the phases at 90°, and 
either keeping all other parameters to their CPX values, or scaling the other parameters 
in such a way that the relation between |mg|, 1^41 and /i is as in the CPX scenario. In 
the following, these ten variants have been studied for rritop = 174.3 GeV/c^ only. 

In all scenarios, theoretical databases provided by the LEP Higgs working group were 
used [7] . In these, radiative corrections have been computed in two different approaches, 
the Feynman-diagrammatic approach of Ref. [30], already selected in the CP-conserving 
case (see Section 4), and the renormalization group approach of Ref. [39]. As in Sec- 
tion 4, the Feynman-diagrammatic calculations use version 2.0 of the FeynHiggs code. 
The renormalization group corrections rely on the CP- violating version CPH of the 
SUBHPOLE code.^ Contrary to the CP-conserving case, where the calculations in 
the Feynman-diagrammatic approach were the most complete due to the inclusion of all 
dominant two-loop order terms, in the case of CP-violation neither of the two calcula- 
tions can be preferred on theoretical grounds. Both contain one and two-loop corrections, 
but the CPH code has a more complete phase dependence at the two-loop order while 
FeynHiggs contains more corrections at the one-loop order with the full complex phase 
dependence and more corrections at the two-loop order but without the full phase de- 
pendence. This may result in large differences when convoluted with the experimental 
inputs. We thus present our results in the two frameworks separately. A comparison 
between the two calculations in the CP-conserving case can be found in Ref. [42]. 

The phenomenology of the three neutral Higgs bosons of the CP- violating MSSM is 
illustrated in Fig. 15 in the case of the CPX scenario for a top quark mass of 174.3 GeV/c^, 
with radiative corrections computed in the renormalization group approach. The top 
figures show that the two lightest scalars. Hi and H2, are likely to be kinematically 
accessible at LEP2 in wide regions of the parameter space, in which their widths remain 
lower than 1 GeV/c^, that is below the experimental resolution. The cross-section curves 
show that at low and large tan/3, the dominant production processes are the HiZ and 
H1H2 processses, respectively, as in the CP-conserving case (see Fig. 4). On the other 
hand, at intermediate tan/5 and moderate mni, the HiZ cross-section is significantly 
weakened, as a result of the suppressed HiZZ couphng due to CP-violation. In the same 
region, the H2Z process compensates only partly for this loss. Finally, the figures showing 
branching fractions compare the dominant Hi and H2 branching fractions for different 
values of tan /?. For all values of tan (3, decays into bb and r+r" saturate the width of the 
lightest Higgs boson. Hi, in the mass range above the bb threshold up to the maximal 
sensitivity of LEP. In the same mass range, the second hghtest Higgs boson, H2, decays 

Since this work, updated versions of the two codes, CPsdperH [40] and FeynHiggs 2.5 [41] have been made available. 
In both cases, the changes concern the Higgs boson decays and have no substantial impact on the phenomenology at LEP. 




JO 



2?, 



00 120 140 



160 180 200 220 240 
mj^2 (GeV/c^) 




60 80 100 120 
mj^j (GeV/c^) 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




tanB 

: 4.6 




CPX scenario 


- 2.9 
: 2.9-- 








4.6"""---.._ 








/; H,^bb 






"' '; ■ '; H,^TT 

;,' '; H^^bb 




4.6 


y ; H2 — >HjH^ 


- 2.9-^ 

4.6 
: 2.9 '" 







10 20 30 40 



50 60 70 80 90 100 
nij^j (GeV/c^) 



- — 

> 
O 

5 



10 



10 



:CPX scenario 






Tjjj 

r 

^ H2 












:20.7- 

10.3—-^;;:^^^^^^ 






38.9 ^^S^^^^"^^""^ 
- 30.2^^^ 




1 


20.7^ 






- 10.3 
_ tanP 







, ii 



34 



20 



40 



60 



80 



100 



120 

2^ 



(GeV/c") 



> 

o 

o 



no 



10 



CPX scenario 


tanp = 20.7 

----- tan|3=10.3 












H|Z 






HjZ \ 




H,Z '\-V 





20 



40 



60 



80 



120 

niui 



100 

'HI 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




: tanP 
: in^ 


CPX scenario 


: 20.7 




L1O.3 




; 10.3 — 






\ H|^bb 




\ Hj^bb 






- 20.7 

-20.7 

: 10.3 









10 20 30 40 50 60 70 80 90 100 

nijjj (GeV/c ) 



Figure 15: Properties of the three neutral Higgs bosons of the CP-violating MSSM in 
the CPX scenario with mtop = 174.3 GeV/c^. Top: Hi, H2 and H3 masses and Hi, H2 
widths at various tan/5 values. Middle: production cross-sections at ^/s = 206 GeV, 
as a function of and tan/5. Bottom: Hi and H2 dominant branching fractions as 
a function of m-Hi and tan (3. Hi and H2 decays into bb (solid and dash-dotted lines) 
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predominantly into bb at large tan/3 only. At low and intermediate tan/?, the cascade 
decay H2 — > HiHi dominates over the bb final state at masses up to 50 GeV/c^ or so. A 
loss in experimental sensitivity can thus be anticipated in regions where the HiZ cross- 
section is negligible and the H2Z signals arc not significant with respect to background, 
due to too weak H2Z cross-sections or H2 branching fractions into fermions. 

6.2 Scan procedure 

The scan procedure is similar to that described in Section 4.2 for the CP-conserving 
scenarios. The only changes are the following. The scan was performed over the MSSM 
parameters tan/3 and mH±. The range in mH± spans from 4 GeV/c^ up to 1 TeV/c^. 
Values of mH± below about 100 GeV/c^ were noticed to give unphysical negative mass 
squared values in most scenarios and thus were removed from the scans. The range in 
tan/3 extends from the minimal value allowed in each scenario up to 40, a value above 
which the Higgs-bottom Yukawa coupling calculation becomes unreliable in the CP- 
violating MSSM scenarios. Theoretical points were generated randomly in both tan/3 
and mH± with a granularity which is sufficient to map the general features of the exclusion 
regions. 

The signal expectations in each channel were then computed as outlined in Section 4.2, 
except for the width effects. In the CP-violating MSSM scans, the widths remain well 
below the experimental resolution for tan/3 up to 40 and below 120 GeV/c^ (see 
Fig. 15). Signal efficiencies and PDFs were thus exclusively determined from simulations 
with Higgs boson widths below 1 GeV/c^. 

7 Results in CP-violating MSSM scenarios 

The regions of the MSSM parameter space excluded at 95% CL or more by combining 
the neutral Higgs boson searches of Table 1 are hereafter discussed in turn for each 
scenario. The additional constraint from the limit on the Z partial width that would 
be due to new physics [35] (used as described in Section 3.2 of [7]) brings no gain in 
sensitivity in any of the scenarios tested. Results are presented only in the (mni, tan/3) 
plane, which is the only one relevant at LEP since the minimal values of and tan /3 
in most scenarios are such that the region accessible at LEP is much reduced in the other 
projections. 

7.1 Dependence on the phases 

The excluded regions in the (mHj, tan/3) plane for the CPX scenario and its vari- 
ants with different phase values are presented in Fig. 16 for the renormalization group 
approach [39] and in Fig. 17 for the Feynman-diagrammatic calculations [30]. The top 
mass value is 174.3 GeV/c^ in all plots. 

Going from 0° to 180°, the excluded domain varies significantly. The quahtative 
trend, valid in the two theoretical approaches, is as follows. The extreme values (0° and 
180°) correspond to scenarios with no CP violation, and hence to a large excluded region. 
Moreover, at 180°, the theoretically allowed region is reduced, especially at large tan /3 due 
to unphysical values of the Higgs-bottom Yukawa couphng. At phases between 60° and 
135°, losses in sensitivity are observed at large tan /3 and above 50 GeV/c^, as well as 
in the intermediate tan (3 range for below 60 GeV/c^. This is the consequence of the 
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Figure 16: CP-violating MSSM scenarios with corrections as in Ref. [39] for different 
values of the phases: regions excluded at 95% CL by combining the results of the neutral 
Higgs boson searches in the whole DELPHI data sample (light-grey). The dashed curves 
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Figure 17: CP-violating MSSM scenarios with corrections as in Ref. [30] for different 
values of the phases: regions excluded at 95% CL by combining the results of the neutral 
Higgs boson searches in the whole DELPHI data sample (light-grey). The dashed curves 
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strong suppression of the HiZZ coupling due to CP-violation, as already encountered in 
Fig. 15 in the case of the CPX scenario. More generally [37], the HiZZ coupling decreases 
slowly (by a few tens of %) with phases below about 75° and is strongly suppressed (by 
three to four orders of magnitude) for phases around 90°. For phases above 100°, the 
coupling is partially restored, mostly at low tan/?. This explains the evolution of the 
upper bound of the experimentally excluded area as a function of phases in Figs. 16 
and 17. The changes are moderate for phases up to 60° and significant for the 90° and 
135° phases, where the experimental sensitivity relies mainly on the H2Z process at low 
tan/? and on the H1H2 production at large tan/?, both giving large signals at moderate 
rriHi only, typically below 60 GeV/c^ (see Fig. 15). 

At the 90° and 135° phases, there are also unexcluded areas at masses lower than 
60 GeV/c^ in the intermediate tan/3 range, between about 4 and 16. These are related 
to weakened sensitivities in the H2Z or H1H2 searches. To take the CPX scenario as 
an example, at masses below 15 GeV/c^, the dominant final state is the (H2 HiHi)Z 
channel. The lack of experimental searches at LEP2 for such final states with be- 
low the bb threshold (see Tab. 1) explains the unexcluded area which is observed at 
these masses, in agreement with the expected sensitivity. The largest value of CLg in 
this region is 52% in the renormalisation group framework and 50% in the Feynman- 
diagrammatic approach. Still in the CPX scenario, the hole at around 50 GeV/c^ 
arises in the region where the decays H2 Hi Hi and H2 — > bb become approximately 
equal, leading to a loss of significance of the H2 signals, as pointed out in Section 6.1 (see 
Fig. 15). The largest value of CLg in this region is 17% (expected 4%) in the renormalisa- 
tion group framework and 37% (expected 11%) in the Feynman-diagrammatic approach. 
In both frameworks, these CLg values are observed at tan/?~4, mHi~50 GeV/c^ and 
mH2~105/107 GeV/c^. The observed exclusion in this region is weaker than expected, 
which is due to a slight excess of data over the expected background. The value of 
1-CLb at the point of weakest exclusion is indeed 15% (corresponding to a 1.4 sigma 
deviation) in the renormalisation group framework and 12% (1.5 sigma deviation) in the 
Feynman-diagrammatic approach. Conversely, the largest deviation in the whole hole 
has a value of 1-CLb of 3.3% (2.1 sigma deviation) in the two approaches. This value 
is observed at tan/?~16, mHi~45 GeV/c^ and mH2~107 GcV/c^ in the renormalisation 
group framework and tan/?~ll, mHi~52 GeV/c^ and mH2'~^lll GeV/c^ in the Feynman- 
diagrammatic approach. At this point, the CLg values are 5.4% (expected 0.1%) and 6.5% 
(expected 0.2%) in the two frameworks, respectively. The combined LEP data show also 
deviations in this region [7]. 

Finally, differences between the two theoretical frameworks are visible mainly at large 
tan/?, where the Feynman-diagrammatic calculations predict significantly higher HiZ 
residual cross-sections (e.g. a factor about 4 in the CPX scenario for between 40 
and 80 GeV/c^), leading to a better experimental sensitivity. Differences in the phase 
dependence of the results are also visible, which are likely to refiect the different phase 
treatment between the two calculations. 



7.2 Dependence on /i and Mgusy 

The excluded regions in the (mnj, tan/?) plane for the CPX scenario and its variants 
with different values of n and Ms^sy are presented in Fig. 18 for the renormalization group 
approach [39] and in Fig. 19 for the Feynman-diagrammatic calculations [30]. In all plots, 
the common CP- violating phase is 90° and the top mass value is 174.3 GeV/c^. 
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Figure 18: CP-violating MSSM scenarios with corrections as in Ref. [39] for different 
values of /i and Msnsy- regions excluded at 95% CL by combining the results of the neutral 
Higgs boson searches in the whole DELPHI data sample (light-grey). The dashed curves 
show the median expected limits. The medium-grey areas are the regions not allowed by 
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Figure 19: CP-violating MSSM scenarios with corrections as in Ref. [30] for different 
values of /i and Msusy: regions excluded at 95% CL by combining the results of the neutral 
Higgs boson searches in the whole DELPHI data sample (light-grey). The dashed curves 
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The dependence of the results on the value of is as expected from the scaling of the 
dominant CP- violating effects with Im{iiA) (see Section 6.1). The exclusion is almost 
entirely restored for values of ji lower than 2 TeV/c^, the value in the CPX scenario, 
and gets weaker at 4 TeV/c^. In the first two variants, despite the CP-violating phase 
being at 90°, there are always two production processes with significant rates in every 
point of the kinematically accessible parameter space. In the variant at 4 TeV/c^, due 
to the large value of ji and the CP- violating phase at 90°, the H1H2 and H2Z processes 
are suppressed for all values of tan/3, as well as the HiZ process at intermediate and 
large tan/3 values. In the Feynman-diagrammatic approach, the HiZ cross-sections are 
partly restored at large tan j3 which explains the difference between the results in the two 
theoretical frameworks in that region. Note also that the theoretically allowed region 
is much reduced at large tan (5 in this scenario due to unphysical values of the bottom 
Yukawa coupling. 

The dependence on the value of Mgusy is presented in the two bottom plots of Figs. 18 
and 19. The first scenario corresponds to setting Mgusy at 1 TeV/c^, twice its value in 
the CPX scenario. As the dominant CP- violating effects are proportional to M~^^^, the 
exclusion is restored in this variant. The reason is as in the case of the two variants 
with low values of i.e. there are always two production processes with significant rates 
in every point of the kinematically accessible parameter space. In the second scenario, 
Msusy is still set at 1 TeV/c^ but the values of |mg|, /x and \A\ are also scaled by a factor 
2, leaving the CP-violating effects almost unchanged (see Section 6.1) with respect to 
the CPX scenario. This explains why the exclusion region in this variant is close to that 
in the CPX scenario. The few differences between the excluded regions in these two 
scenarios are due to different cross-sections for some of the processes which contribute 
most to the experimental sensitivity, that is the H1H2 and H2Z processes at masses below 
60 GeV/c^, and the HiZ process at higher masses. As an example, the better coverage 
of the low mass region at intermediate tan fi values in the scaled variant is explained by 
slightly higher H1H2 cross-sections. 

7.3 Dependence on mtop 

The excluded regions in the (mn^, tan/3) plane for the CPX scenario with different 
mtop values are presented in Fig. 20 for the renormalization group approach and in Fig. 21 
for the Feynman-diagrammatic calculations. 

The results show a strong dependence on the value of mtop, as expected since the 
dominant CP- violating effects scale with m^^p. In the two theoretical approaches, the 
exclusion in the intermediate tan/3 range is gradually reduced as mtop increases and 
eventually vanishes for tan/3 between about 3 and 5 and a top mass of 183 GeV/c^. 
This can be traced to the suppression of the H2Z and H1H2 cross-sections with increasing 
values of mtop, leaving no significant rate in any of the three possible production channels. 
At large tan/3, as mtop increases, the H1H2 cross-section is reduced and the exclusion 
gets weaker in the renormalization group approach while it is almost unchanged in the 
Feynman-diagrammatic framework. As already mentioned, this is a consequence of the 
higher HiZ residual cross-sections predicted by the latter calculations at large tan/3. 

7.4 Summary 

Scans of the CPX scenario and its variants revealed that CP violation in the Higgs 
sector can have a significant impact on the experimental sensitivity of LEP. The strong 
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Figure 20: CP-violating MSSM scenarios with corrections as in Ref. [39] for different 
values of mtop: regions excluded at 95% CL by combining the results of the neutral Higgs 
boson searches in the whole DELPHI data sample (light-grey). The dashed curves show 
the median expected limits. The medium-grey areas are the regions not allowed by theory. 
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Figure 21: CP-violating MSSM scenarios with corrections as in Ref. [30] for different 
values of mtop: regions excluded at 95% CL by combining the results of the neutral Higgs 
boson searches in the whole DELPHI data sample (light-grey). The dashed curves show 
the median expected limits. The medium-grey areas are the regions not allowed by theory. 
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suppression of the neutral Higgs boson couplings to the Z boson translates into a loss of 
redundancy in the different search channels, and hence leads to a reduced coverage of the 
parameter space. The most significant reduction is observed in the intermediate tan/3 
region, typically between 3 and 10, down to the lowest Hi masses. It occurs for phases 
between 90° and 135°, top mass values equal to 174.3 GeV/c^ or larger, and values of 
the ratio \iiA\/ M^^^y equal to 8 or larger. As a consequence, no absolute mass limits 
can be derived in these scenarios. On the other hand, the low tan /3 region appears still 
disfavoured, as in the CP-conserving models. Scans were performed using two different 
theoretical approaches for the radiative correction calculations. Although the two sets of 
results show large differences, they both lead to the same qualitative conclusions. 

8 Conclusions 

Searches for Higgs bosons in the whole data sample of the DELPHI experiment have 
been combined to derive constraints on MSSM benchmark scenarios, including models 
with CP-violation in the Higgs sector. Experimental results encompass searches for 
neutral Higgs bosons in dominant final states expected in most MSSM models, as well 
as searches for charged Higgs bosons and for neutral Higgs bosons decaying into hadrons 
of any flavour, which bring a gain in sensitivity in restricted regions of the parameter 
space. An additional improvement is obtained by applying the experimental results to 
more production processes than the two expected main channels, namely the associated 
production of the lightest Higgs boson with a Z boson and the pair-production of the 
two hghtest Higgs bosons. In the CP-conserving MSSM, the experimental sensitivity 
at LEP relies on the hZ, hA and HZ channels, the last leading to a significant gain in 
sensitivity in scenarios where the third neutral Higgs boson, H, is kinematically accessible. 
In the CP-violating MSSM, the total signal at LEP is spread mainly over the HiZ, H2Z 
and H1H2 channels. Accounting for the simultaneous production of all possible signals is 
essential in this type of scenario where CP-violating effects can lead to strong suppression 
of one channel or another. 

In all CP-conserving scenarios, the experimental results allow a large fraction of the 
parameter space to be excluded, even in scenarios designed to test potentially difficult 
cases (e.g. vanishing production cross-sections or decay branching fractions) either at 
LEP or at hadron colliders. Limits on masses of the h and A bosons were deduced as 
well as upper and lower exclusion bounds in tan/3. The dependence of these limits on 
mtop was studied in a range between 169.2 to 183.0 GeV/c^. To quote but one result, 
the following hmits at 95% of CL have been established in the framework of the m™^ 
scenario with mtop — 174.3 GeV/c^: 

mh> 89.7 GeV/c^ and mA> 90.4 GeV/c^ for any tan/3 between 0.4 and 50, 
tan/3<0.72 or tan/3>1.96 for any mA between 0.02 and 1000 GeV/c^. 

These mass limits are insensitive to variations of the top quark mass. The excluded 
range in tan /3 decreases with increasing mtop and would vanish if mtop was as large as 
183.0 GeV/c^. This scenario provides the most conservative bounds on tan/3 among the 
eight CP-conserving scenarios tested. 

In the CP-violating scenarios, large domains of the kinematically accessible parame- 
ter space remain unexcluded due to strong suppressions of the couplings between the Z 
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and the Higgs bosons induced by CP-violation. Hence no absolute limits can be set on 
the Higgs boson masses in these scenarios. The unexcluded areas arise in the intermedi- 
ate tan (5 range, typically between 3 and 10. Their contours vary considerably with the 
value of mtop and the MSSM parameters which govern the CP-violating effects, 
Msusy and the phase arg{A). These scenarios have been studied in two different theoret- 
ical frameworks for the radiative correction calculations. The impact of CP-violation is 
observed to be qualitatively the same in the two approaches. 
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Appendix 1 

We give hereafter efficiencies of the Yukawa T+T~bb analysis pubhshed in Ref. [21] 
and apphed here to the hA -^r'^r'qq signal. 
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Table 7: hA — >T+T~qq channel : efficiencies of the selection (in %) at LEPl as a function 
of the masses of the A and h bosons. The analysis, described in Ref. [21], was designed to 
search for Yukawa production in the T"'"T~bb final state. The quoted errors are statistical 
only. 
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Appendix 2 

We give hereafter efficiencies of the h — > qq analyses pubhshed in Ref. [20,1] and 
apphed to (h — > AA— >cccc)(Z — > qq) signals with low A masses. 
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Table 8: (h — > AA)(Z — > qq) channel with A — > cc: efficiencies of the selection (in %) at 
v^s = 199.6 and 206.5 GeV as a function of the masses of the A and h bosons, for mA 
between the cc and bb thresholds. Efficiencies at higher masses can be found in Ref. [20,1]. 
We refer the reader to Ref. [1] for the definition of the two operational periods of the 
2000 data taking campaign. The quoted errors are statistical only. 
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